Laser cooling all the way down to Fermi superfluid 
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Numerical simulations for realistic experimental parame- 
ters demonstrate that laser cooling on the attractive side of 
the Feshbach resonance can drive fermions much below the 
superfluid transition. For the assumed set of experimental 
parameters the transition takes place at 0.352>, and laser 
cooling can drive the system down to at least 0.0857> in a 
time of a few seconds. Superfluid growth is self-consistently 
included in simulations. 
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Fermi superfluids are at the forefront of both experi- 
mental and theoretical activity in the physics of ultracold 
dilute atomic gases. In addition to being a challenge in 
their own right, they provide promising model systems to 
test the theory of fermionic superfluidity, or the Bardeen- 
Cooper-Schrieffer (BCS) theory, in the regime of strong 
interactions, when the size of Cooper pairs becomes com- 
parable, or even less than the average interparticle dis- 
tance. In this regime excitations of the system develop a 
pseudogap indicating that Cooper pairs are quasi bound 
states of pairs of fermions [1]. This is also the limit that 
seems to be a key to understanding the physics of high T c 
superconductors, where pseudogap effects have been ob- 
served experimentally. The cold atoms provide a model 
system of high flexibility thanks to application of the Fes- 
hbach resonance technique [2]. The Feshbach resonance 
between an open and closed channel allows one to manip- 
ulate an effective s-wave scattering length between pairs 
of atoms with opposite "spin" . This technique allows 
one to drive the system through the BCS-BEC crossover, 
when the scattering length diverges and changes sign 
from negative (attractive) to positive (repulsive) . At the 
same time the state of the system continuously evolves 
from the BCS superfluid of weakly correlated Cooper 
pairs, through the pseudogap regime, and then the uni- 
tarity limit, when the scattering length diverges, to the 
Bose-Einstein condesate of biatomic molecules. 

The main problem in the experiments on atomic Fermi 
superfluids is how to cool the system sufficiently deep into 
the quantum degenerate regime to reach the superfluid- 
ity. Unfortunately the standard evaporative cooling stops 
to work in the degenerate regime, but before the superflu- 
idity is achieved [3] . One of the reasons may be creation 
of holes deep in the Fermi sea by collisions with the buffer 
gas [4]. However, this problem can be circumvented in an 
ingenious way thanks, again, to the Feshbach technique 



[5]. The idea is to obtain first a molecular condensate 
on the BEC side of the Feshbach resonance, and then 
continuously transform it into BCS superfluid by an adi- 
abatic passage to the other side of the resonance. This 
adiabatic transition reduces the temperature sufficiently 
to reach the superfluid state. By now, several succesful 
experiments with molecular condesates were performed, 
and paved a way to first succesful experiments with Fermi 
superfluids [6,7]. The Innsbruck experiment [7] is the first 
to report direct signatures of Fermi superfluidity in quasi- 
particle excitation spectrum, in excellent agreement with 
the theory of Refs. [8]. 

The purpose of this Letter is twofold. First, we demon- 
strate that it is possible to laser cool a Fermi gas down 
to the superfluid BCS state on the attractive (BCS) side 
of the Feshbach resonance. This is challenging because, 
so far, the standard laser cooling has not even allowed 
for observing quantum degeneracy. However, theoretical 
studies of laser cooling show in particular that it is pos- 
sible to reach the degeneracy in Fermi systems [9] . Laser 
cooling offers an alternative route towards BCS super- 
fluid directly on the negative side of the Feshbach reso- 
nance. Moreover, it opens a possibility to cool the BCS 
superfluid prepared via the "standard" method down to 
even lower temperatures. The second goal of this Let- 
ter is to develop for the first time fully self-consistent 
quantum kinetic theory of BCS state creation via laser 
cooling . The results are very promising: despite the 
fact that both BCS formation and normal mean field ef- 
fects shift the energies of quasiparticlc excitations with 
respect to the equally spaced harmonic oscillator levels, 
the cooling remains efficient at least down to the temper- 
ature T ~ 0.085Ti?. This temperature is well below the 
superfluid transition temperature which for the assumed 
reasonable experimental parameters is close to 0.35Tp. 
Optimization of the cooling protocol could significantly 
reduce this temperature. 

We consider two species of fermions with "spin up 
and down" , molecules which are their bound states, and 
fermions in an excited state. Populations of the spin up 
and down fermions are the same. We assume that mu- 
tual s-wave interaction between these two populations 
leads to thermalisation on a timescale, which is faster 
than the rate of laser cooling. The system remains in 



1 similar to self-consitent theories of the BEC growth [10] 
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a quasi-cquilibrium state with slowly decreasing temper- 
ature. This equlibrium state is described by the gen- 
eralised weak coupling BCS theory, or boson-fermion 
model [11] which includes Bosc-Einstcin condensate of 
the molecules. Structure of the ground state and its ex- 
citations follows from the set of Bogoliubov-de Genncs 
equations together with self-consistent definitions of the 
gap function and the mean field potential. Laser cooling 
of fermions was described in detail in Ref. [9]. Coher- 
ent laser excites atoms from, say, the spin down ground 
state to the excited state and at the same time sponta- 
neous emission brings them back to the spin down ground 
state. Frequencies of the laser assure that the generalized 
Raman cooling takes place [12]. In this Letter we take 
into account interactions between fermions and describe 
the cooling process in terms of instantaneous Bogoli- 
ubov quasiparticlcs, whose eigenfunctions and eigenen- 
ergies are self-consistently updated during the evolution. 
We work in the festina lente limit to avoid reabsorbtion 
effects [13], and also employ spherical symmetry and er- 
godic approximations (see also [10]). 

The boson-fermion model [11] is defined by the Hamil- 
tonian 



H= Jdh 



+ 0o$-$lM++ 

=+,-,e 

+ h.c. + {2v - 2n)fiJI> m . (1) 



The fields ip+ and -0_ describe the fermions with spin up 
and down respectively, ip e corresponds to the fermions 
in the excited state, and ip m is the bosonic molecular 
field. Here Hi = — |^V 2 + V(r) — p, is a single particle 
Hamiltonian, m is atomic mass, V(r) is a harmonic trap 

potential, p is chemical potential, 50 = 47r ^ a " is bare 
interaction strength with a negative bare s-wave scatter- 
ing length fl , A is a coupling between pairs of atoms 
and molecules, and v is detuning from the Feshbach res- 
onance. 

In the mean field approximation, which closely fol- 
lows the BCS theory, the quartic interaction term in 
the Hamiltonian (1) is made quadratic by replacing 
all products of pairs of operators with their averages. 
The averages are the mean field potential W{r) = 

go (tp±(r)ip±(r) ), and the anomalous potential of the 



pairing field P{r) = yp + (r)tp-(r)j , which is mixing 
fermions with spin up and down. The molecular field 
4> m is replaced by the amplitude of molecular condensate 
(<j) m ). After ellimination of the molecular condensate 
one obtains linear equations for the spin up and down 
fermions 



Hi $± + W(r) 4,± T g P{?) , (2) 



with an effective coupling constant g 



go 



These equations mix tfj and 1/^, but can be "diagonalizcd" 
by the Bogoliubov transformation 



i>+(r) = b„ h+ u rn {r) - P m _ v* m {r) , 



V>_(r) 



(r) + 



(3) 
(4) 



with fermionic quasiparticle annihilation operators 6 m ,±. 
The Bogoliubov modes (u m ,v m ) fulfill the Bogoliubov-dc 
Gennes equations 



Hi u r 
Ho v„ 



W u, 

W V„. 



A iv 



A* Um , 



(5) 



with positive energies u m of quasiparticle excitations. 
Here A = — gP is the quasiparticle gap function. In 
a thermal state with inverse temperature (3, the aver- 
age occupation numbers of quasiparticle states are given 
by Fermi-Dirac distribution N m — (S^ ^fem^) = (1 + 
exp(/3w m )) _1 . Equations (5) are solved together with 
the self-consistency conditions 



A(f) 



-gJ2(l-™ m )u m (r)v* m (r) 



(6) 
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W(f) = ff ]T(l - N m )\v m (r)\ 2 + N m \u m (r)\ 2 , (7) 

m 

by succesive iterations with the chemical potential p ad- 
justed to keep the average number of atoms constant. 
The ultraviolet divergence in Eq.(6) is regularized using 
the quickly convergent method of Ref. [14]. 

Excitation of atoms from the spin down state to the 
excited state is described by the Hamiltonian 



las 



-Qe ikLF ipl(r)ip^{r) + h.c. - 5^ e 



(8) 



driving coherent oscillations with Rabi frequency O and 
laser detuning 8. The excitation is accompanied by the 
spontaneous emission, which after rotating wave approx- 
imation is described by a superoperator 

£ P = T X! u ™kD[b* m _e k ]p + V mk T>[b m ,+ek}p, (9) 



mk 



Here the 



with a spontaneous emission rate of 7. 
Lindblad superoperator is T>[A]p = ApA* — 
^A'Ap, and the matrix elements are e.g. U m k — 
J drtkW(flk)\u m k(k)\ 2 with the spontaneous emission 
pattern W(flk) and the generalized Frank-Condon fac- 
tors u m k{k) = J d 3 r e lkr wl(r)u m (f) . Here Wk{r) is the 
fc-th eigenstate of the harmonic oscillator and an an- 
nihilation operator of an excited atom in this state. 
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Adiabatic cllimination of the excited state, similar as in 
Ref. [15], leads to the kinetic equations for the occupation 
numbers taking into account the effects of Fermi statistics 



dN, 



(U = E C 1 - N ^~) N n,- -(n~m) + 

n 

C mn (1 - JV m ,_) (1 - AT„ ;+ ) - A nm N m ^N n!+ 



dN, 



dt 



= E r ™-« (! - N ^+) N ",+ - ( n " m ) + 

n 

C„ m (1 - N m , + ) (1 - JV n ,_) - A m „JV m>+ AT ni _ . 



They describe relaxation of ± quasiparticles (r^i_ n ), and 
creation/annihilation of pairs of + and — quasiparticles 
(C nm and A nm ). The transition rates are e.g. 



It") = 



27 



E 



7 2 ^mfc|M»fc(fcL)| 2 

(<5-^+u;„) 2 +7 fe 2 



(10) 



Here tu^ is the energy of the fc-th harmonic oscillator 
state and 7fc is approximate spontaneous decay rate of an 
excited atom in this state: jk = lJ2 m U m k{l — N m -) + 
V m kN mj+ . 

Laser cooling drives average occupation numbers N m ,± 
of quasiparticle states out of thermal equilibrium. At 
the same time interactions between quasiparticle states 
drive the system towards thermal equilibrium. Numeri- 
cal simulations in Ref. [9] show that for cooling rates in 
ms range thermal relaxation remains faster than cool- 
ing all the way down to T ~ 0.03 T F . This jus- 
tifies our assumption of fast equilibration to a quasi- 
cquilibrium state. After each short period of laser 
cooling the occupation numbers N mi ± go out of equi- 
librium, where N m = [1 + cxp(/3cj m )] _1 . The ini- 
tial total energy E(T) = ± u m N mt ± changes by 
dE = Em^rn {N m ,+ + N m> l - 2N m ). The relax- 
ation brings the occupation numbers N mj ± to a new 
state of equilibrium at a temperature T + dT, but it does 
not change the total energy. The energy of the system 
at the new temperature E(T + dT) differs only slightly 
from the initial E(T) by dE T = E(T + dT) - E(T) w 
Em w i^m(l — N m ). Conservation of energy in ther- 
mal relaxation means that dE = dEx- In our simulations 
we use this equality to find the new lower temperature 
T + dT, and then solve the Bogoliubov-de Gennes equa- 
tions self-consistently to adjust the Bogoliubov modes 
and the chemical potential to the lower temperature. 
With the new Bogoliubov modes we calculate new tran- 
sition rates etc. 

Calculation of the transition rates is the most time 
consuming part of the numerical simulation. This is why 
we were forced to assume spherical symmetry. For spher- 
ically symmetric A(r) and W(r) the Bogoliubov modes 
can be separated as e.g. w/„(r) Y lm (9, <p) and their ener- 
gies tuin do not depend on m. To be consistent with the 



spherical symmetry we also made the ergodic approxima- 
tion [9,10], assuming that quasiparticle occupation num- 
bers Ni n ^± do not depent on m. In other words, we as- 
sume fast equilibration within each quasiparticle energy 
shell, which keeps the system in a spherically symmetric 
state. 

In our simulations we use 81 harmonic oscillator lev- 
els (shells) and assume that there are N = 10660 atoms 
with spin up and down. In the non-interacting case at 
T = the atoms fill 39 energy shells. The frequency of 
the isotropic trap is w — 2ir 2400Hz. The bare scatter- 
ing length for the interaction between the two species is 
ao = — 157<zb, where cib is the Bohr radius. This value 
corresponds to the interactions between the spin states 
\F = 9/2, m F = 9/2} and \F = 9/2, m F = 7/2) of 40 K. In 
the harmonic oscillator units this scattering length gives 
go = —0.32. In our calculations we use two values for 
the effective interaction strength: g = —0.32 and g = — 1 
enhanced by the Feshbach resonance. The wavelength of 
the cooling laser is A = 720nm and the laser detuning 
is S = —12Ulo. In all simulations the Rabi frequency is 
much less than the spontaneous emission rate, Q = O.I7, 
so that average occupation of excited state remains small. 
This is one of the assumptions in the adiabatic cllimina- 
tion of the excited state. 

In Fig.l we show temperature as a function of time for 
laser cooling with four values of 7 consistent with the fes- 
tina lente regime. Each plot in this figure consists in fact 
of two indistinguishable plots, one for the bare g = —0.32 
and the other for the enhanced g = — 1. For the stronger 
g = — 1 there is a transition to the superfluid state at low 
temperature (see Fig. 2), but the non-zero gap function 
which appears at the transition does not influence the 
efficiency of laser cooling. The process of laser cooling is 
dominated by the higly degenerate states with large an- 
gular momenta I w Ep/fiui whose quasiparticle energies 
are not influenced much by the gap function localized in 
the center of the trap [8,16]. In more physical terms, the 
superfluid droplet in the center of the trap is surrounded 
by a shell of normal Fermi gas. As long as this meso- 
scopic effect persists the laser cooling in the superfluid 
state remains as efficient as in the normal state. 

The initial cooling is the fastest for 7 = 10 but the 
lowest temperature of 0.0857V is obtained for 7 = 1.25. 
The best strategy to reach low temperatures in reason- 
able time is to do the cooling in a few stages starting with 
large 7 and finishing the job with small 7, or to continu- 
ously decrease 7 as the temperature decreases. The inset 
in Fig.l shows an example of a four stage cooling process. 

Summarizing, we made numerical simulations for rea- 
sonable experimental parameters consistent with the re- 
quirements of the festina lente regime which demonstrate 
that it is possible to laser cool fermions on the nega- 
tive side of the Feshbach resonance all the way down 
to the Fermi superfluid. Despite the pairing effects and 
the usual mean field interaction, which shift quasiparti- 
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cle energies with respect to the equally spaced harmonic 
oscillator levels, it is experimentally risible to reach at 
least 0.085Tp in a time as short as a few seconds. 
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FIG. 1. Temperature as a function of time during the laser 
cooling for four values of 7: 7 = 10 (solid lines), 7 = 5 
(dashed lines), 7 = 2.5 (dashed-dotted lines), and 7 = 1.25 
(dotted lines). Each pair of plots in this figure corresponds 
to g = —0.32 (higher temperature) and g = — 1 (lower tem- 
perature). The lowest temperatures of 0.079Tf and 0.085Tf 
respectively can be reached for 7 = 1.25. On the other 
hand, the initial stage of laser cooling is efficient at most for 
7 = 10, when the temperature of 0.15Tf is reached within 
0.85s. The inset shows laser cooling to the superfluid state 
with g = — l realized in four stages to reach the lowest tem- 
perature in reasonable time. The spontaneous emission rate 
is switched as 7 = 10 — > 5 — > 2.5 — > 1.25 at the temperatures 
T/T F = 0.16,0.11,0.10, correspondingly. 
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FIG. 2. The maximal magnitude of the gap function in the 
center of the trap \A(r — Q)\/Thu for the same simulations 
with g — —1 as in Fig. 1. The inset shows the gap in the four 
stage cooling of the inset of Fig.l. 
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